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ABSTRACT 

We present the detection of 34 Lya emission-line galaxy candidates in a 80 x 80 x 60 co- moving Mpc 
region surrounding the known z — 2.38 galaxy cluster J2143-4423. The space density of Lya emitters 
is comparable to that found by Steidel et al. when targeting a cluster at redshift 3.09, but is a factor 
of 5.8 ± 2.5 greater than that found by field samples at similar redshifts. 

The distribution of these galaxy candidates contains several 5-10 Mpc scale voids. We compare our 
observations with mock catalogs derived from the VIRGO consortium ACDM n-body simulations. 
Fewer than 1% of the mock catalogues contain voids as large as we observe. Our observations thus 
tentatively suggest that the galaxy distribution at redshift 2.38 contains larger voids than predicted 
by current models. 

Three of the candidate galaxies and one previously discovered galaxy have the large luminosities and 
extended morphologies of "Lya blobs" . 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: fundamental parameters 



1. INTRODUCTION 

The two-point correlation coefficient of high redshift 
(z > 2) galaxies is quite similar to that of galaxies today 
(e.g. Steidel et al. 1996, 2000, Giavalisco et al. 1998). 
The topology of the distribution of high redshift galax- 
ies, however, is not yet clear. The two-point correlation 
coefficient alone is not very sensitive to this topology. 
In the local universe, a large fraction of galaxies lie in 
filaments and sheets, such as the Great Wall (Geller & 
Huchra 1989). These filaments may be several hundred 
Mpc in length, and are separated by voids which can 
be tens of Mpc across. Our understanding of how and 
when filaments and voids were established depends on 
the geometry of the universe and biasing of galaxies (e.g. 
White & Frenk 1991; Kauffmann, Nusser, & Steinmetz 
1997; Mo, Mao, & White 1999). 

Mapping such large scale structure at high redshifts 
requires a careful balance of field of view (to encompass 
the structures), redshift coverage (to accept enough test 
objects in the structure while avoiding confusion by over- 
lapping structures in distance) , and photometric sensitiv- 
ity (to find enough test objects). Photometric redshifts 
(see Hogg et al. 1998 for a review) provide an efficient es- 
timate of redshift for large imaging surveys, but they are 
not accurate enough to avoid confusion from multiple 
overlapping structures, and cannot delineate the struc- 
tures directly. Bright QSOs may be markers of dense re- 
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gions (Ellingson, Yee & Green, 1991), but are too sparse 
to provide multiple samples within individual structure 
elements. 

Narrow-band imaging to find Lya emitting galaxies has 
many advantages as a technique for mapping the topol- 
ogy of the galaxy distribution at high redshifts. The tech- 
nique picks out galaxies in a narrow-range of redshifts, so 
the two-dimensional distribution of their positions on the 
sky can be used to constrain the topology alone, without 
the need for expensive follow-up spectroscopy At suffi- 
ciently low flux limits the space density of detectable Lya 
emitters becomes comparable to that of Lyman-break 
galaxies (eg. Hu & McMahon 1996). They should be 
dense tracers of large scale structure. With large-format 
prime-focus imagers, it is relatively straightforward to 
survey very large volumes of the early universe. The fast 
beam required for large field coverage puts a lower limit 
the width of the narrow-band filter. Typically this lim- 
its the line-of-sight depth covered in an exposure to be 
roughly equal to or greater than the width. A disad- 
vantage of Lya searches is that the Lya emission from 
a galaxy is very hard to predict theoretically, due to the 
very high optical depth in this line, its dependency on 
metalicity and star formation rate, and the ease with 
which it can be obscured by dust. 

Two studies find some evidence for filamentary struc- 
ture in the distribution of Lya emitting galaxies on scales 
of around 5Mpc (Campos et al. 1999, M0ller & Fynbo 
2001). Other surveys have found evidence for clustering 
in Lya emitting galaxies, but not for filamentary struc- 
ture (Steidel et al. 2000, Ouchi et al. 2002). 

In this paper, we present the detection of 34 candidate 
Lya emitting galaxies in a region containing J2143-4423, 
a z=2.38 cluster of galaxies and damped Lya absorbers 
(Francis et al., 1996, 1997, 2000). Our survey covers a 
region 80x80x60 co- moving Mpc in size, and hence is 
sensitive to much larger structures than previous stud- 
ies at comparable redshifts. The narrow-band imaging 
also detects four extended nebulae of Lya emission, of- 
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Table 1. Observations 



Filter 


cxptime / frame 
(s) 


no. frames 


FWHM 

(") 


5a 1 
(AB mag.) 


U 


1800 


5 


1.49 


23.9 


B 


600 


12 


1.37 


26.2 


V 


600 


12 


1.37 


25.3 


R 


600 


5 


1.64 


24.0 


I 


600 


5 


1.11 


23.8 


NB 


1800,2700 


5.1 


1.35 


23.5 



The 5<t limit for a point source, on the AB magnitude system. 



ten dubbed "blobs" (Steidel et al., 1999, Keel et al. 1999, 
Francis et al., 2000). Throughout the paper, we assume 
a A-dominated flat Universe (Ho — 65 km s" 1 Mpc -1 , 
Q. M = 0.3, fi A = 0.7). 

2. OBSERVATIONS AND DATA REDUCTION 

The data were taken with the MOSAIC II instrument 
on the Blanco 4m telescope at CTIO. Observations were 
made on the nights of 7-8 August, 1999, in the Johnson 
UBV and Cousins RI filters, and in a custom narrow- 
band interference filter (NB) centered at 4107A with a 
full width at half maximum (FWHM) of 54A. The NB fil- 
ter was designed to image Lya emission at z=2.38 (Fran- 
cis et al. 2000). The delivered filter covers Lya redshifts 
between 2.356 < z < 2.401. The Mosaic II is a prime- 
focus camera consisting of 8 CCDs, each with 2048x4096 
pixels. The camera has a plate scale of 0.27"/pixel and 
covers an overall area of (36') 2 with small 9-13" gaps 
between the CCDs. Individual exposures were dithered 
by 10'. The full NB integration time of 19,800 sec was 
achieved for a 26'x26' region centered on the J2143-4423 
cluster, while shorter integrations extend over a 48'x50' 
field. The final trimmed images are 45.9'x45.5' in size 
excluding 9'x9' in the South-East corner with only one 
exposure. The observations are summarized in Tabled 

This was the first science run with the MOSAIC II. 
The readout clocking was still being tuned and biases 
had a slowly varying harmonic pattern due to the beat- 
ing of unsynchronized clocks. An FFT based algorithm 
was developed to subtract the ringing component of the 
bias. The period of the pattern was approximately 160 
pixels (120") in the readout direction on the CCD. Be- 
cause this scale is significantly larger than the size of any 
source in the images the subtraction algorithm should 
not affect the detection or photometry of of the sources. 
All detections were verified in the raw images. 

Further image reductions were performed using the MS- 
CRED package in IRAF 9 . These include a correction for 
crosstalk between the CCDs, bias subtraction and flat- 
fielding. 

Twilight flats were used for all of the images except for 
the I-band. The I-band twilight flats exhibited consid- 
erable fringing with a different pattern than that due to 
the night sky. The I-band dome flats however were not 
uniformly illuminated and an illumination correction was 
derived from highly smoothed twilight flats. The I-band 
images were defringed by combining all of the images 

9 IRAF is distributed by NOAO, which is operated by AURA 
Inc., under contract to the NSF. 



with object rejection. The objects were rejected by mak- 
ing a mask using the SExTRACTOR (Bertin & Arnouts 
1996) object images. The object mask images were con- 
volved with a 20 pixel diameter top hat kernel to reject 
the low surface brightness tails of the objects. 

Wide-field optical distortions were corrected using a 
standard astrometric solution for MOSAIC II provided 
in IRAF. Final alignment of the images was performed 
by matching each image to approximately 2000 stars in 
the APM catalog. The alignment was good to 0.3 pixels 
(0.08") rms. 

Before combining the images we created bad pixel 
masks which included hot pixels bad columns and cos- 
mic rays. Cosmic rays were selected using the jcrrejl 
package in IRAF (Rhoads 2000). We found this pack- 
age to work extremely well for finding individual cosmic 
rays, however we had difficulty balancing the parameters 
to eliminate the faint halos associated with the cosmic 
ray hits. To eliminate these we devised a filtering al- 
gorithm which calculated the Tookey biweight statistic 
(Beers, Flynn & Gebhardt 1990) in a region around each 
CR and eliminated all points adjacent to the CR which 
were > 2.5cr from the estimated background level. 

Photometric calibration was performed using 56 stars 
in the Landolt fields 90, 92, SA110 and 113 (Landolt 
1992). A small correction was include in the calibration 
to account for Galactic extinction (Schlegel, Finkbeiner 
& Davis 1998). Catalogs of individual sources were com- 
piled using SExTRACTOR. The limiting magnitudes vary 
across the images. Figure Q shows an exposure map for 
the NB imaging. The map was created by summing a set 
of the flat field images scaled by exposure times and with 
the same offsets as the NB exposures. The 5cr detection 
limits for point-like objects in each color are summarized 
in Table The NB filter was used as a primary de- 
tection band, with measurements made in other bands 
using the same apertures, in order to find objects with 
weak continua but strong lines. Areas that are covered 
by only one exposure, such as the south-east corner are 
excluded from the analysis. Photometry was measured 
in the Kron aperture (Kron 1980), defined as 2 times the 
first moment of the radial light distribution. The first 
moment is approximately equal to the half light radius 
of the distribution. The photometric error was measured 
by SExtractor. Table 1 lists the 5cr photometric sensitiv- 
ity for a point source in each band. 

3. RESULTS 

We detect 2450 point-like objects and 7000 resolved 
objects above the 5ct limit in the NB filter. Lya emitters 
can be identified in this sample by excess flux in the NB 
filter compared to the continuum B filter. The limits for 
NB — B color excess were determined by modeling the 
distribution of photometric errors in the NB — B color 
as a function of the NB magnitude (see Teplitz et al. 
1998). The outer envelope of 5er errors that defines the 
detection limit for emission-line candidates is fit by the 
function: 

[{NB — 15.25)/7.6] 4 + 0.15 (1) 

However, bright emission-line candidates with low 
equivalent width (EW) are more likely to be low red- 
shift interlopers than cluster members (see section 13.10 . 
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Fig. 1. — Exposure map for the NB imaging. North is up East 
is to the left. 



Fig. 2. — Number Counts for resolved objects detected in the 
narrow-band images. Assuming the number counts can be extrapo- 
lated by a power law the detections are 50% complete at NB=23.25 



We consider candidates above the EW = 125 A limit to 
be likely Lya emitters. This limit is obtained as a broad- 
minus narrow-band color through the relation: 



NB -B = -2.5 log 



1 + EW/Wb 
1 + EW/Wj 



NB 



(2) 



where Wb and Wmb are the widths of the broad- and 
narrow-band filters. 

In addition, objects that are only marginally detected 
in the narrow-band are unlikely to be good candidates. 
For simplicity, we establish a uniform narrow-band mag- 
nitude cut rather than base the cut on the varying depth 
across the image. The magnitude cut is based on the 
50% completeness limit for narrow-band detections (see 
Figure EJ). 

Figure |3 shows the NB — B color excess for all objects 
detected in the NB image. The Lya emitters from Fran- 
cis et al. (1996) are clearly detected in the new NB mea- 
surement. For new emission line candidates we consider 
objects that are detected above the 50% completeness 
limit in the narrow-band filter (NB < 23.5), are located 
above the 5a color excess line, and have an observed EW 
greater than 125 A. We can detect galaxies with emis- 
sion line strengths brighter than 1.4 x 10~ 16 ergs cm~ 2 
s . This flux limit corresponds to a Lya luminosity 
of 1.9 x 10 9 L Q . Thirty-seven spatially resolved NB ex- 
cess objects are detected. Tabled lists the positions and 
measured properties of these objects. 

In addition, 7 unresolved objects with U BV colors con- 
sistent with QSO's (e.g., Hall et al. 1996a, 1996b) are 
detected in emission with equivalent widths greater than 
30A in the observed frame. Table 01 lists the positions 
and measured properties of these objects. 

Finally, a QSO at z=2.38 was found by Hawkins 
(2000). We detect this object in emission, but it was 
not selected as a QSO candidate because the UBV col- 
ors are within the stellar locus. We include this QSO in 
Table 




Fig. 3. — The excess flux in the NB filter as a function of mag- 
nitude. Black points are extended sources and grey are unresolved 
sources. The large filled circles are the three spectroscopically con- 
firmed Lya-emitters (Francis et al. 1997). Upward pointing arrows 
are objects undetected in the B-band filter. Three damped Lya 
absorbers (Francis & Hewett 1993, Francis, Wilson & Woodgate 
2001) are indicated by open stars. The curved lines indicated the 
5<r confidence level for excess emission; that is, points above the 
line are candidate emission-line object. The horizontal line is the 
EW b s = 125 A level. The vertical line is the 50% completeness 
cutoff; that is, only candidates brighter than that line are consid- 
ered likely to be real. Objects which are simultaneously above the 
5(7 threshold, above the observed EW b s = 125 A level and brighter 
than the 50% completeness limit are selected as Lya-emission can- 
didates. 



3.1. Foreground Contamination 

For spatially resolved narrow-band excess sources the 
only likely contaminating line is [O II] (3727A). Our 
narrow-band filter centered at 4107A, is sensitive to [O II] 
emission from foreground 0.095 < z < 0.109 galaxies. In- 
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Table 2. Emission-Line Galaxies 



a (2000) 


<5 (2000) 


NB 


B - NB 


Lya flux 


EW 


Notes 






(mag.) 


(mag.) 


(erg cm -2 scc _1 ) 


(A) 




21:40:19.90 


a a in /ion 

-44:19:48.2 


22.62 


1.91 


n TOwin — I f i 

2.73x10 ' 


380 




ni ,i n oi nn 

21:4U:ol.UU 


a a o/" 1 n a cr 

-44:oo:U4.5 


22.97 


1.53 


i oiwin — IB 

1.81x10 


215 




21:40:33.10 


/i /i o -in o 

-44:36:10.8 


22.91 


1.37 


i o i v / 1 n — 1 fi 

1.81x10 


169 




ni /in oc 'vt 

21:40:30.77 


/i A nn no o 

-44:2U:2o.o 


23.20 


1.78 


1 r/ 1 i / 1 n — 1 fi 

1.56x10 


312 




m /in ao nn 

21:40:48.09 


A A O "1 n "1 T 

-44:31:01.7 


23.11 


2.87 


i n r v / 1 n — I fi 

1.95x10 


3026 




m a n jo nT 

21:40:48.97 


a A n "i no c 

-44:01:23.6 


22.77 


1.29 


n n n , , 1 n — I fi 

2.00x10 


150 




m /in r" o nn 

21:40:58.22 


a a nn nn n 

-44:00:22.0 


23.00 


3.00 


n 1 o v / 1 n — If! 

2.18x10 ° 


5900 




m j "i nn nn 

21:41:02.90 


/i /i n "i p" r n 

-44:01:55.9 


22.81 


1.25 


i on,, i n — I fi 

1.89x10 


141 




ni j "i n t oo 

21:41:07.38 


/I /I O O "1 "1 T 

-44:38:11.7 


23.04 


1.34 


i p" n , , 1 n — I f i 

1.59x10 


162 




m /ii /i/i a i 

21:41:44.41 


/i /I T7 nr t 

-44:37:06.7 


22.68 


1.82 


o tT/iwin — ItS 

2.54X1U u 


331 




21:41:47.66 


/i /i n "i n "i n 

-44:21:21.9 


22.61 


2.37 


n n p" , , 1 n — If) 

2.95x10 


824 




21:41:5o.55 


/I /I O O 1 o o 

-44:38:18.3 


23.18 


1.82 


i £?nwin — 16 
l.OUXlU 


331 




n 1 i o n/ -1 no 

21:42:Uo.Uo 


A A O A A T C\ 

-44:o4:47.9 


23.18 


2.58 


i Tn w 1 n — 1 fi 

1.79X1U u 


1277 


B6 


n 1 /in i/i no 

21:42:14.28 


a a on "i r" o 

-44:32:15.8 


22.18 


2.07 


/i m , , i n — In 

4.21x10 


489 




n 1 /in nT c 

21:42:27. 5o 


a a nn on 1 

-44:2U:oU.l 


20.97 


2.47 


i otTwin — IS 

l.o5xl(J 


1004 


Bl 


m /in no C/i 

21:42:28. 54 


/i /i on o o cr 

-44:32:38.5 


23.20 


2.44 


i TOwin — Ifi 

1.7oXl(J u 


944 




n 1 /in nn to 

21:42:29. 7o 


,1 /i ni nn o 

-44:21:U2.8 


22.86 


2.30 


n o o w 1 n — Ifi 

2.ooXl(J 


723 


B2 


n 1 i o on nn 

21:42:o2.2U 


/i a nn "i o i" 

-44:2U:lo.o 


22.90 


2.53 


n on,, 1 n — Ifi 

2.oUXlU 


1141 


B4 


m /ino/i oo 

21:42:o4.8o 


a A nT n/" 1 n 

-44:27 :Uo. 2 


21.80 


3.00 


nr,, in — Ifi 

D.95X1U 


5900 


B7 


n 1 /in /in/^o 

21:42:42.63 


/i /i o n n n n 

-44:30:09.0 


21.20 


3.00 


i 1 ,i,,in — "IS 

1.14x10 


5900 




n 1 /in c ^ nT 

21:42:54. U7 


a a 1/1 on T 

-44:14:o9.7 


22.93 


1.37 


i TOwin — Ifi 

1.7oXl(J u 


169 




n 1 /in r"/^ 0,1 

21:42:56.34 


/i /i o t r" o 

-44:37:56.8 


22.48 


1.56 


n o z' 1 , , 1 n — If) 

2.86x10 


225 




n 1 /i o nn nn 

21:43:00.09 


/i /I 1 n ni T 

-44:19:21.7 


22.57 


1.70 


o to w 1 n — Ifi 

2.7oXlU 


277 




m /to no C7 

21:43:03.57 


/i /i no /i /i n 

-44:2o:44.2 


21.42 


2.61 


n n/""win — Ifi 

9.UbXlU u 


1371 


B5 


n "i /i o no on 

21:43:03.80 


/I /I O 1 A A f\ 

-44:31:44.9 


22.16 


2.30 


A A O v , 1 n — If) 

4.43x10 


723 




n 1 /tone nn 

21:43:05.90 


a a nT ni n 

-44:27:21. U 


21.06 


2.10 


i 1 n w 1 n — IS 
1.19X1U 


513 




n 1 /ion/^ /in 

21:43:06.42 


a a nT nn c 

-44:27:00.6 


22.48 


2.61 


o /ii, , 1 n — Ifi 

3.41x10 


1371 




21:43:11.48 


-43:59:01.0 


23.13 


2.87 


1.91xl0~ 16 


3026 




21:43:22.22 


-44:13:06.5 


23.09 


1.93 


1.78X10- 16 


392 




21:43:23.80 


-44:41:36.4 


22.87 


1.79 


2.12X10" 16 


317 




21:43:24.06 


-44:27:59.9 


22.42 


1.81 


3.22xl0~ 16 


326 




21:43:37.41 


-44:17:53.4 


23.21 


2.51 


1.73xl0~ 16 


1092 




21:43:37.48 


-44:23:52.8 


22.65 


2.47 


2.88xl0~ lfi 


1004 




21:43:44.92 


-44:05:46.3 


22.53 


1.66 


2.81xl0~ 16 


261 




21:43:48.30 


-44:08:26.9 


22.42 


1.48 


2.95xl0~ 16 


200 




21:44:12.15 


-44:05:46.6 


21.46 


2.41 


8.56xl0~ 16 


890 




21:44:12.97 


-43:57:56.7 


22.70 


1.53 


2.31X10- 16 


215 





Table 3. Emission-Line QSO Candidates 



a (2000) <5 (2000) NB B - NB line flux EWq B-V U - B Notes 

(mag.) (mag.) (erg cm _2 scc _1 ) (A) (mag.) (mag.) 



21:40:52.46 


-44:36:21.35 


19.96 


0.83 


2.04e-15 


70 


0.26 


-0.28 


21:41:54.50 


-44:18:35.97 


20.14 


1.01 


1.96c-15 


96 


0.05 


-0.31 


21:42:35.12 


-44:32:30.36 


21.35 


0.99 


6.35e-16 


93 


0.36 


-0.11 


21:42:43.49 


-44:14:25.32 


20.14 


0.41 


1.02c-15 


27 


0.07 


-0.35 


21:43:16.44 


-44:16:51.54 


21.04 


0.67 


6.51c-16 


51 


0.08 


-0.30 


21:43:20.38 


-44:20:00.35 


20.59 


0.46 


7.39e-16 


31 


0.06 


-0.52 


21:43:22.53 


-44:31:49.18 


19.88 


0.98 


2.45e-15 


91 


0.04 


-0.43 


21:43:26.25 


-44:26:03.44 


20.77 


0.49 


6.58e-16 


34 


0.26 


-0.33 



deed, we detect emission from many well resolved spiral 
galaxies below our EW cutoff. 

[0 11] interlopers within our Lya sample would have 
absolute B magnitudes around —13.8, and would hence 
be Blue Compact/HII galaxies (BCGs). At this redshift 
they would only be marginally spatially resolved, so wc 
could not separate them from z = 2.38 Lya emitters by 
morphology alone. 

To be detected through our filter, they must lie within 
a box of comoving size 5.9x5.9x60 Mpc, and have a rest- 



frame [O II] equivalent width greater than 114A. Using 
the local luminosity function of Jerjen, Bingelli & Free- 
man (2000) we would expect to find ~ 12 BCG galaxies 
in such a box. Pustilnik ct al. (1999), however, found 
that only ~ 8% of BCG galaxies have [O II] equivalent 
widths large enough to meet our selection criteria. If this 
foreground region was an average one, we would therefore 
expect ~ 1 foreground BCG galaxy to be contaminating 
our sample. 

The southernmost 80% of our field was included in the 
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Fig. 4. — The distribution of emission-line objects in the field. 
The small open circles represent emission-line galaxy candidates. 
The large filled circles are extended emission-line nebulae (blobs). 
The large open stars are the three known DLAs, and the small 
open stars are quasar candidates at the cluster redshift or lower 
redshifts. East is to the left and North is upwards. 



Las Campanas Redshift Survey (LCRS, Shectman et al. 
1996), which has excellent sensitivity to large scale struc- 
tures at redshift z ~ 0.1. The region in which foreground 
[O II] emitters could contaminate our survey contains 
a single isolated cluster (Abell 3800, Abell, Corwin & 
Olowin 1989), but is otherwise empty of galaxies. The 
LCRS detects many filaments of galaxies at this redshift, 
but none lie in our field at redshifts that could cause 
[O II] emission to impersonate Lya at z=2.38. Despite 
the presence of Abell 3800, the galaxy density is indeed 
slightly below the LCRS average at this redshift. Abell 
3800 is located 10' west and 3' south of our field cen- 
ter (see Fig . We see no concentration of candidate 
Lya emitting galaxies at this location: indeed quite the 
opposite. Only one of our candidates lies within one pro- 
jected Mpc of this location. Three others lie around 1.3 
projected Mpc from the cluster, but these are the three 
for which we have spectroscopic confirmation that they 
lie at z = 2.38. We therefore tentatively conclude that 
foreground contamination is unlikely to be a big problem. 

Foreground contamination is more likely a problem for 
the 7 narrow-band excess QSO candidates. The narrow- 
band excess could be caused by Lya emission at redshift 
2.38, but could also be quasars at other redshifts, with 
different lines producing the excess. Indeed one of these 
sources was found (Francis et al. 1997) to be a QSO at 
z=1.66, with its C IV (1549A) emission producing the 
narrow-band excess. The luminosity and redshift dis- 
tribution of QSOs implies that about 60% or 4 of these 
quasars should be Lya emitters at redshift 2.38 (see Palu- 
nas et al. 2001). If some of these sources are QSOs at 
redshift 2.38, it will be further evidence that QSOs lie in 
the most massive dark matter halos, which in turn exist 
in the most highly clustered environments (e.g. Silk & 
Weinberg 1991). 

4. DISCUSSION 



RA (arcmin) 

Fig. 5. — Density contours for the distribution of Lya candidates 
estimated with the adaptive kernel routine. The lowest contour is 
one galaxy per resolution element 3.3 X 10 — 3 galaxies arcmin -2 . 
The highest contour is 5.7 X 10 -2 galaxies arcmin -2 . 



4.1. The Space Density of Lya Emitters 

We are sensitive to Lya emitting galaxies over a 
45'x45' field minus 9'x9'in the South-East corner, in 
the redshift range 2.356 < z < 2.401. For our adopted 
cosmology, this corresponds to a volume 6100 co-moving 
Mpc 2 in area and 59 co-moving Mpc deep. Figure \S\ 
shows the density contours of Lya candidates estimated 
using an adaptive kernel routine (Beers, Flynn & Geb- 
hardt 1990). The initial smoothing scale is 12 arcmin- 
utes. The highest density contour is 5.7 x 10~ 2 sources 
per square arcminute corresponding to a spatial density 
of 3.0 x 10~ 4 sources per cubic Mpc. The average den- 
sity over the whole field is 1.9 x 10 -2 sources per square 
arcminute corresponding to a spatial density of 1 x 10~ 4 
sources per cubic Mpc. 

Steidel et al. (2000) surveyed 20,000 cubic co-moving 
Mpc at redshift 3.09, targeting a cluster of Ly-dropout 
galaxies. Their survey was sensitive to lower Lya lu- 
minosities than ours, and they used a lower rest-frame 
equivalent width threshold. They find 9 candidates meet- 
ing our luminosity and equivalent width thresholds cor- 
rected for redshift, yielding a density of 4.5 x 10~ 4 sources 
per cubic Mpc. Steidel et al. calculate that this space 
density is a factor of 6.0±2.4 greater than that of smaller 
survey of Cowie & Hu (1998). This implies that peak 
density in our field is a factor of 4 overdense. 

4.2. Star Formation Rates 

To estimate the star formation rate we first assume 
that star formation powers the observed Lya flux. This 
is almost certainly a great underestimate: even tiny 
amounts of dust will greatly reduce the observed Lya 
flux due to its high optical depth. We use an un-reddened 
Lya/Ha ratio of 8:1, and the conversion factor between 
Ha flux and star formation rate of Kennicutt (1983). 
Our faintest candidates (NB=23.25) have inferred star 
formation rates of ~ 9M0yr -1 . 

Summing over all galaxies down to our flux limit, 
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and making a crude correction for incompleteness 
amongst the faintest, the integrated star forma- 
tion rate density (assuming no dust obscuration) is 
2.4 x 10~ 3 M Q Mpc~ 3 yr- 1 : well below the value (~ 
O.lM0Mpc~ 3 yr -1 ) derived from integrated blue light by 
Madau, Pozzetti & Dickinson (1998). 

We also estimate star formation rate from the 1600A 
UV continuum using the Madau et al. formula. At 
z=2.38 the 1600A continuum falls in the U-band. We 
detect 23 of our candidates in the V-band. The star- 
formation density implied by these candidates is 5 x 
10~ 3 A/ Q Mpc~ 3 yr~ 1 . The star formation rate for a can- 
didate at the V-band limit of 25.3mag is ~ 18M yr _1 . 
The 14 galaxies not detected in the V-band could add, 
as an upper limit, 0.6 x lO _3 M0Mpc~ 3 yr _1 to the star 
formation rate density. 

Our Lya candidates therefore contribute only a small 
fraction (~ 5%) of the overall star formation rate density 
at z=2.38. 

In a detailed study of the spectra of Lyman-break 
galaxies, Shapley et al. (2003) conclude that galaxies 
with the strongest Lya emission have bluer UV continua, 
weaker interstellar absorption and smaller star-formation 
rates than galaxies with weak Lya emission or Lya ab- 
sorption. Lya emission from galaxies with the highest 
star-formation rates is absorbed by dust. 

4.3. The Distribution o/Lya Emission Candidates 

The distribution of our candidate resolved Lya emit- 
ting sources is shown in Fig^J There is no obvious over- 
density surrounding the previously known cluster at the 
field center. Instead, most of the candidates lie in a broad 
band extending from the north-east of the field to the 
south. There is a significant is the lack of galaxies in 
the north and west. These regions have long integration 
times. There are three close pairs or triplets of candi- 
dates, but otherwise the candidates are well spaced. The 
pairs have a spacing of about 20 arcsec or a minimum sep- 
aration of 200 proper kpc. The spacing of the galaxies in 
the triplet is about 50 arcsec or a minimum separation 
of 440 proper kpc. 

Comparison to a random distribution 

We used two statistics to compare our data to a ran- 
dom distribution. We compare these statistics to results 
from Monte Carlo simulations generated with the same 
exposure-time mask. In this section we show that the 
distribution of our candidates has a significant excess of 
close (< 1') pairs, and of large (6' - 8') voids. 

The first statistic is essentially the angular two-point 
correlation function. For a series of angular separations, 
we take each data point in turn and count the number 
of other data points within that angular radius. 

The second statistic is a two-dimensional analogue of 
the void probability function (VPF). For each angular 
scale, we randomly place 1000 circles with that radius 
on our field, so that the circles do not extend beyond the 
edge of our data. We count the fraction of these circles 
which did not contain any data points. The requirement 
that the circles not extend beyond the edge of the data 
makes this statistic most sensitive to voids near the cen- 
ter of the field. Void probability functions are sensitive 
to structures such as voids and filaments, to which the 
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Fig. 6. — The void probability function of our data (solid tri- 
angles) compared to the average void probability functions of ran- 
domly distributed points (solid line) and void probability functions 
derived from the ACDM n-body simulations (open squares). 



two-point correlation function is relatively insensitive. 

Our data were compared to random Monte Carlo sim- 
ulations generated as follows. At each location in our im- 
age, the total exposure time was calculated. Using these 
total exposure times, the relative flux limits in every part 
of our field were computed. The number/magnitude dis- 
tribution of our brighter candidates can be reasonably 
approximated as a power-law, with the number per unit 
magnitude increasing by a factor of roughly 5. Combin- 
ing this number-magnitude relation with the flux lim- 
its, we calculated the relative probability of finding a 
candidate per unit area in any given part of our image. 
We then used this relative probability map to generate 
2000 fake data sets, each containing the same number of 
sources as in the real data. 

Our data have a significant excess of galaxies with sep- 
arations less than 1'. Only 0.1% of the Monte Carlo sim- 
ulations produced this many close galaxies. On larger 
scales, our data still show an excess of galaxies, but this 
is no longer significant at the 5% level. 

The void probability function (Fig EJ) is more inter- 
esting. There is a clear excess of voids on scales of 5 
arcminutes or greater. On scales of 5' - 8', this excess is 
significant at the 99.9% confidence level (ie. fewer than 
0.1% of our randomly generated data sets had VPFs as 
large as our data). 

We conclude that there is a significant excess of voids, 
and of small scale clustering, in our data. This conclu- 
sion is insensitive to the details of how we generated the 
random data sets: changes in (for example) the assumed 
number/magnitude relation made little difference to our 
results. 

Comparison with XCDM n-body Simulations 

The distribution of our Lya candidates is thus inconsis- 
tent with a random distribution, with 99.5% confidence. 
Assuming that they are indeed galaxies at redshift 2.38, 
is their distribution consistent with theoretical predic- 
tions? In this section, we show that while the two-point 
correlation coefficient of our galaxy candidates is quite 
consistent with the predictions of a ACDM simulation, 
the distribution of our candidates shows more voids on 
6' - 8' scales, with 97% confidence. 

To evaluate the apparent excess of voids, we compared 
our measured distribution of galaxies against the ACDM 
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n-body simulations performed by the VIRGO consortium 
(Kauffmann et al. 1999). We used their galaxy catalog, 
generated for redshift 2.12. For each galaxy, they list its 
position, stellar mass, gaseous mass and star formation 
rate. 

Mock catalogs were generated as follows. Our observa- 
tions were targeted at the source Bl, a probable pair of 
giant elliptical galaxies surrounded by an extensive neu- 
tral hydrogen cloud (Francis et al. 2000). Bl has an 
estimated stellar mass of 1O 1O ' 9 M0 and gaseous mass of 
10 9 - 7 M Q . We picked five galaxies in the ACDM simu- 
lations that had stellar and gaseous masses at least this 
large, and which lay far enough from the edge of the sim- 
ulation volume. The regions around these "Bl-analogue" 
galaxies should be good matches to our data. For com- 
parison, we also picked eight random locations within the 
ACDM simulation volume. We then extracted the galax- 
ies that lay within a volume 70x70x46 co-moving Mpc 
(for our adopted cosmology) centered at each of these 
locations. There were 180 ± 30 (1 er) galaxies in their 
catalogs in each of these volumes. 

We then had to sparsely sample the galaxies within 
each of these volumes to generate mock catalogs with 
(a) 37 galaxies, and (b) a probability of being detected 
consistent with our exposure-time map. We did this in 
three different ways. 

1. Randomly. As Lya emission can be generated 
in many ways (by star formation, shocks, cool- 
ing flows, AGN etc), and its escape in measurable 
quantities depends on accidents of the geometry 
and dust distribution inside galaxies, we first as- 
sumed that all the galaxies in the ACDM simula- 
tions had the same probability (37/180 ~ 20%) of 
being detected by our survey. 

2. Star Formation Rate. If the Lya flux is generated 
by star formation, only galaxies with star forma- 
tion rates in excess of 8Af yr _1 would have been 
detected. We therefore took only galaxies with star 
formation rates at least this large, and generated 
random sub-samples with the correct size and prob- 
ability of detection as a function of position in our 
field. Only 25% of the galaxies within the ACDM 
cubes had star formation rates this large. 

3. Mass. We took only galaxies with masses greater 
than 10 10 ' 3 Mq, once again picking random sub- 
samples of the appropriate size and with proba- 
bilities of detection scaling appropriately with ex- 
posure time in each part of our image. 

For each data cube, and for each of these three 
sub-sampling techniques, we repeated the random sub- 
sampling 2000 times. We then calculated the angular 
two-point correlation function and the void probability 
function for each sub-sampling of each data cube, and 
compared our results with the observations. 

Fig[7|shows our data, and five of the mock ACDM data 
cubes (centered on the five Bl-analogue galaxies). The 
cubes shown have been randomly sub-sampled: those 
sampled on the basis of star formation rate or mass look 
very similar. Fig shows the average void probability 
function of the randomly sub-sampled data cubes (cen- 
tered on Bl-analogue galaxies). 
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Fig. 7. — The distribution on the sky or our candidate 
Lyaemitting galaxies (top left) compared to five simulated data 
sets drawn from the ACDM n-body simulations. These simula- 
tions use random sub-sampling and are centered on Bl-analogue 
galaxies. 



The angular two-point correlation coefficient of our 
data, and of the ACDM simulations are in excellent agree- 
ment. As Fig shows, however, the agreement is not as 
good for the void probability function. Our data have a 
void probability function considerably higher than that 
of the simulations on large scales. 

Is this difference significant? The VPFs of our 13 sim- 
ulated data cubes (five centered on Bl-analogue galax- 
ies and eight randomly located) actually showed rather 
little scatter, when averaged over all the 2000 random 
sub-samplings of each. The dominant source of scatter 
was the selection of the 37 galaxies in each simulated vol- 
ume, which was typically twice as large. Our uncertainty 
is thus dominated by the small sample size rather than 
by cosmic variance. 

If we take the five data cubes centered on Bl- 
analogues, and generate 2000 random sub-samplings of 
each, fewer that 0.1% of all these random sub-samplings 
have void probability functions as large as we observe on 
scales of 5' - 8'. If we sub-sample on the basis of star for- 
mation rate or galaxy mass, this fraction rises to around 
1%. We would expect this, as more massive galaxies, and 
those with larger star formation rates, are expected to be 
more strongly biased. If we repeat these calculations for 
the randomly centered data cubes, the void probability 
function also rises slightly. This is presumably because 
centering our data cube on a Bl analogue precludes the 
possibility of a void in the center of the field. 

Thus the observed distribution of galaxies has an ex- 
cess of voids over the ACDM simulations, significant at 
roughly the 99% level. A qualitative impression of this 
can be gained from Fig [7J While all the ACDM simula- 
tions show empty regions near the edge of the field, this 
is an artifact of the shorter exposure times in these re- 
gions. The simulations do not show large voids in the 
central ±13 arcminutes. Due to the restriction that we 
only place voids so that they do not overlap the edges 
of our field, our VPF statistic is mostly sensitive to the 
distribution of galaxies in this central region. 

Discussion of the Voids 

We thus have tentative evidence from our survey of 
larger voids in the galaxy distribution at redshift 2.38 
than predicted by one particular ACDM simulation. 



With only 97% confidence, a sample size of 34 and no 
spectral confirmation for most of our candidate galaxies, 
this evidence must be regarded as tentative at best. But 
if it is confirmed by larger samples, what is it telling us? 

Steidel et al. (2000) did not see voids, but they cov- 
ered too small a region on the sky to be sensitive to 
the scales of void we are detecting. Ouchi et al. (2003) 
imaged a region ~ 25% the size of ours, searching for 
Lya candidates at z = 4.86. Voids can be seen in their 
data if only brightest objects are selected so that their 
number per unit area matches ours. Voids might also 
explain the enormous field-to-field variance in the space 
density of Lya sources found by Pascarelle, Windhorst 
& Keel (1998). The spiky distribution of the redshifts 
of Lyman-dropout galaxies seen by Steidel et al. (1998) 
and Adelberger et al. (1998) would also be consistent 
with large voids in the high redshift galaxy distribution. 

If the distribution of galaxies really contains more 10 
co-moving Mpc scale voids than predicted, either ACDM 
is predicting too uniform a distribution of galaxies, or the 
voids contain dark matter halos which for some reason 
do not contain detectable galaxies. 

There is some evidence for the latter hypothesis. The 
sight-line to background QSO 2138—4427 passes through 
one of the voids in our data, but the QSO spectrum 
shows a metal enriched damped Lya absorption-system 
at our redshift (Francis & Hewett 1993, Francis, Wilson 
& Woodgate 2001). So there is clearly at least one dense 
concentration of gas and some stars within one of our 
voids. 

4.4. Blobs 

Some high redshift galaxy clusters contain extended 
Lya emission nebulae, or "blobs" (Steidel et al. 1999). 
The blobs are radio quiet (< 140 /xJy at 1.4 GHz), 
but otherwise share morphological similarity with the 
nebulae around radio galaxies (e.g. Kurk et al. 2001). 
Lya blobs are large, bright gas clouds (~ 100/i — 1 kpc, 
n m = 0.3,fi A = 0.7; Lya flux up to 1.8x ~ 10~ 15 erg 
cm~ 2 s _1 ). They appear to be common in regions of 
significant galaxy overdensity (~ 10 times that of the 
field) at high redshift. Most blobs are associated with a 
host galaxy, though not symmetrically centered on it. 
Blobs may break into smaller knots of Lya and con- 
tinuum emission, with velocity differences Ay < 2000 
km s _1 (Steidel et al. 2000). Blobs can be among the 
brightest high-z sub- mm sources (Chapman et al. 2001). 
Some blobs (Francis et al. 2001) show CIV emission (and 
may host AGN) but others do not. 

We detect four Lya blobs in the narrow-band image. 
One of these, "Bl" , has previously been detected (Francis 
et al. 2001). Figure (SJ shows the broad- and narrow- 
band images of each blob, together with the continuum 
subtracted emission-line image. Each blob conforms to 
the expected size and brightness of this class of object. 
Table 0] summarizes the characteristics of the blobs. 

These blobs may be pre-galactic clouds of gas, in the 
process of initial collapse or destined to be stripped away 
into the intracluster medium (ICM). In the blobs, we 
may witness a stage of structure formation which is illu- 
minated by an unusual or transitory phenomenon. The 
excitation mechanism for the Lya emission is not known, 
but three primary models have been suggested: 

In the early stages of a galaxy's formation, gas col- 
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Fig. 8.— REDUCED RESOLUTION. Images of the blobs. The 
top row shows the broad band image of each blob; the middle row 
the narrow-band image; and the bottom row the difference of the 
two, scaled to give an image of the Lya emission. Each image is 
30" X 30". 



lapses into the dark matter potential well and cools via 
radiation (e.g. Fabian et al. 1986). Haiman, Spaans & 
Quataert (2000) calculate that cooling flows (with gas 
temperatures of a few xlO 4 K) could produce the ob- 
served flux range, size, surface brightness profiles and 
number density observed from Lya blobs. While Fardal 
et al. (2000) doubt that cooling flows would have suffi- 
cient luminosity, Francis et al. (2001) suggest that shocks 
within the infalling gas can trigger the blob emission in 
the cooling flow. Confirmation that the blobs are high 
redshift cooling flows would have profound implications 
for the interpretation of semianalytic galaxy formation 
models (e.g. Fall & Efstathiou 1980) since we would be 
observing the initial collapse of gas into dark matter po- 
tential wells. 

Taniguchi et al. (2001 and references therein) suggest 
that starburst-powered winds trigger Lya emission in 
the expelled gas. Dust could obscure the host galaxy, 
with the Lya emission arising from a superwind extend- 
ing ~ 100 kpc from the host. The dusty host galaxies 
would eventually evolve into ordinary elliptical galaxies, 
and the blobs would presumably dissipate into the intr- 
acluster medium. The SCUBA detections support this 
model, implying rapid star formation (Chapman et al. 
2001) and a probability of hosting an AGN (Barger et al. 
2001). Superwinds would indicate that the blobs are the 
next phase of galaxy evolution after the initial collapse. 
Taniguchi & Shioya (2000) estimate this phase would last 
only 0.1 Gyr, occurring after the first 0.5 Gyr of star for- 
mation. The expelled gas would eventually be stripped 
into the ICM, though it would be metal-enriched having 
been processed by supernovae in the host galaxy. 

An AGN may photo-ionize the gas, similar to radio 
galaxy nebulae and Seyfert 2 galaxies with extended 
emission line regions (e.g. Villar-Martin, Tadhunter & 
Clark 1997). For example, a Lya nebula has been found 
around the HDF-S QSO which is radio-quiet (Bergeron 
et al. 1999; Palunas et al. 2000a). For the blobs, we do 
not see a strong UV continuum source; however, there 
may be an AGN surrounded by an obscuring torus (e.g. 
Krolik & Begelman 1986). In fact, most of the Chan- 
dra point sources do not show optical evidence for AGN 
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Table 4. Lya Blobs 



Name 


Rcdshift 


major axis 


Lya flux 


CIV emission 


Ref. a 




w 


(") 


(erg cm sec -1 ) 






SB-1 


3.09 


17 


1.4 x 10~ 15 


no 


1,2 


SB-2 


3.09 


15 


1.2 x 10~ 15 


no 


1,2 


KB-A 


2.39 


~8 


1.1 x 10~ 15 


yes 


3,4 


KB-B 


2.39 


~8 


0.45 x 10~ 15 


yes 


3,4 


Bl 


2.38 


17 


1.8 x 10~ 15 


yes 


5, this work 


B5 


2.38 


8 


1.5 x 10~ 15 




this work 


B6 


2.38 


8 


1.4 x 10~ 15 




this work 


B7 


2.38 


7 


0.7 x 10~ 15 




this work 



a l - Steidel et al. 2000; 2 - Bautz & Garmirc 2001; 3 - Windhorst et al. 1998; 4 - Keel ct al. 1999; 5 - Francis ct al. 2001 



(Mushotzky et al. 2000). Assuming an ionization pa- 
rameter of 10 and a density of 1 cm -3 at a distance of 
10 kpc, a central source with luminosity 9 x 10 45 ergs 
s _1 is required to ionize a Lya blob with L x ~ 10 44 ergs 
s _1 . The gas in AGN-supportcd blobs would contribute 
to the ICM, but with potentially lower metallicity than 
that from a starburst, thus diluting rather than enriching 
the ICM. 



some environmental dependence on whether a galaxy 
emits detectable Lya emission. 

The detection of too many voids in our data is sugges- 
tive, but further observations with much more uniform 
sensitivity over much wider fields, combined with follow- 
up spectroscopy, will be required to definitively measure 
the topology of the distribution of high redshift galaxies. 



5. SUMMARY 

Narrow-band imaging in the J2143-4423 region has 
revealed 34 candidate Lya emitters, including three 
new extended Lya blob candidates and five possible 
quasars. Together with the previously reported emission- 
line galaxies and DLA systems, these detections suggest 
this field is one of the most highly evolved structures 
detected at high redshift. 

The distribution of the candidates is non-random. 
There are more closely grouped galaxies, and more voids 
in the distribution, than would be expected from either a 
random distribution of galaxies or from the ACDM sim- 
ulations of Kauffmann et al. (1999). This may indicate 



We thank K. Weaver and R. Mushotzky for useful dis- 
cussions, and the VIRGO consortium for making the out- 
puts of their n-body simulations publicly available. We 
thank the referee W. Keel for comments which improved 
the manuscript. This study was funded by a NASA grant 
NRA-98-03-UVG-011. CTIO is operated by the Associ- 
ation of Universities for Research in Astronomy (AURA), 
Inc., under cooperative agreement with the National Sci- 
ence Foundation. 

This work was supported by the STIS IDT through the 
National Optical Astronomical Observatories and by the 
Goddard Space Flight Center. 



REFERENCES 



Abell, G. O., Corwin, H. C, & Olowin, R. P. 1989, ApJS, 70, 1 
Adelberger, K. L., Steidel, C. C, Giavalisco, M., Dickinson, M., 

Pettini, M., & Kellogg, M. 1998, ApJ, 505, 18 
Barger, A. J., Cowie, L. L., Mushotzky, R. F., & Richards, E. A. 

2001, AJ, 121, 662 
Bautz, M. W., & Garmire, G. P. 2001, AAS, 199, 100.10 
Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 
Bergeron, J., Petitjcan, P., Cristiani, S., Arnouts, S., Bresolin, F., 

& Fasano, G. 1999, A&A, 343, L40 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Campos, A., Yahil, A., Windhorst, R. A., Richards, E. A., 

Pascarelle, S., Impey, C, & Petry, C. 1999, ApJL, 511, L4 
Chapman, S. C, Lewis, G. F., Scott, D., Richards, E., Borys, C, 

Steidel, C. C, Adelberger, K. L., & Shapley, A. E. 2001, ApJ, 

548, L17 

Cowie, L. L., & Hu, E. M. 1998, AJ, 115, 1319 

Ellingson, E., Yee, H. K. C, & Green, R. F. 1991, ApJ, 371, 49 

Fabian, A. C, Arnaud, K. A., Nulsen, P. E. J., & Mushotzky, R. F. 

1986, ApJ, 305, 9 
Fall, S. M., & Efstathiou, G. 1980, MNRAS, 193, 189 
Fardal, M. A., Katz, N., Gardner, J. P., Hernquist, L., Weinberg, 

D. H., & Dave, R. 2001, ApJ, 562, 605 
Francis, P. J., & Hewett, P. C. 1993, AJ, 105, 1633 
Francis, P. J., et al. 1996, ApJ, 457, 490 

Francis, P. J., Woodgatc, B. E., & Danks, A. C. 1997, ApJ, 482, 
L25 

Francis, P. J., Wilson, G. M., & Woodgate, B. E. 2001, PASA, 18, 
64 

Francis, P. J., et al. 2001, ApJ, 554, 1001 

Geller, M. J., & Huchra, J. P. 1989, Science, 246, 897 



Giavalisco, M., Steidel, C. C, Adelberger, K. L., Dickinson, M. E., 

Pettini, M., & Kellogg, M. 1998, ApJ, 503, 543 
Haiman, Z., Spaans, M., & Quataert, E. 2000, ApJ, 537, L5 
Hall, P. B., Osmer, P. S., Green, R. F., Porter, A. C, & Warren, 

S. J. 1996a, ApJ, 462, 614 
Hall, P. B., Osmer, P. S., Green, R. F., Porter, A. C, & Warren, 

S. J. 1996b, ApJ, 471, 1073 
Hawkins, M. R. S. 2000, A&AS, 143, 465 
Hogg, D. W. et al. 1998, AJ, 115, 1418 
Hu, E. M., & McMahon, R. G. 1996, Nature, 382, 281 
Jerjen, H., Bingelli, B., & Freeman, K. C. 2000, AJ, 119, 593 
Kauffmann, G., Nusser, A., & Steinmetz, M. 1997, MNRAS, 286, 

795 

Kauffmann, G., Colberg, J. M., Diaferio, A., & White, S. D. M. 

1999, MNRAS ,303, 188 
Keel, W. C, Cohen, S. H., Windhorst, R. A., & Waddington, I. 

1999, AJ, 118, 2547 
Kennicutt, R. C, Jr. 1983, ApJ, 272, 54 
Krolik, J. H., & Begelman, M. C. 1986, ApJ, 308, L55 
Kron, R. G. 1980, ApJS, 43, 305 

Kurk, J. D., Pentericci, L., R0ttgering, H. J. A., & Miley, G. K. 

2001, ApSSS, 277, 543 
Landolt, A. U. 1992, AJ, 104, 340 

Le Fevre, O., Deltorn, J. N., Crampton, D., & Dickinson, M. 1996, 
ApJ., 471, Lll 

Lowcnthal, J. D., Hogan, C. J., Green, R. F., Caulet, A., Woodgatc, 

B. E., Brown, L., & Foltz, C. B. 1991, ApJ, 377, L73 
Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 
Mo, H. J., Mau, S., & White, S. D. M. 1999, MNRAS, 304, 175 
M0llcr, P., & Fynbo, J. U. 2001, A& A, 372, L57 



10 



Mushotzky, R. F., Cowie, L. L., Barger, A. J., & Arnaud, K. A. 

2000, Nature, 404, 459 
Ouchi, M. et al. 2003, ApJ, 582, 60O 

Pascarelle, S. M., Windhorst, R. A., & Keel, W. C. 1998, AJ, 116, 
2659 

Pascarelle, S. M., Windhorst, R. A., Keel, W. C, & Odewahn, S. C. 

1996, Nature, 383, 45 
Palunas, P., et al. 2000, ApJ, 541, 61 
Pritchet, C. J. 1994, PASP, 106, 1052 

Pustilnik, S. A., Engcls, D., Ugryumov, A. V., Lipovetsky, V. A., 
Hagen, H.-J., Kniazev, A. Y., Izotov, Y. I., & Richtcr, G. 1999, 
A & AS, 137, 299 

Rhoads, J. E. 2000, PASP, 112, 703 

Shapley, A. E., Steidel, C. C, Pettini, M., & Adelberger, K. L. 

2003, ApJ, 588, 65 
Schlcgcl, D. J., Finkbeincr, D. P., & Davis, M. 1998, ApJ, 500, 525 



Shcctman, S. A., Landy, S. D., Oemler, A., Tucker, D. L., Huan, 

L., Kirshner, R. P., & Schechter, P. L. 1996, ApJ, 470, 172 
Silk, J. & Weinberg, D. 1991, Nature, 350, 272 
Steidel, C. C, Giavalisco, M., Pettini, M., Dickinson, M., & 

Adelberger, K. L. 1996, ApJ Letters 462, L17 
Steidel, C. C, Adelberger, K. L., Dickinson, M., Giavalisco, M., 

Pettini, M., & Kellogg, M. 1998, ApJ, 492, 428 
Steidel, C C, Adelberger, K. L., Shapley, A. E., Pettini, M., 

Dickinson, M., & Giavalisco, M. 2000, ApJ, 532, 170 
Taniguchi, Y., & Shioya, Y. 2000, ApJ, 532, L13 
Taniguchi, Y., Shioya, Y., & Kakazu, Y. 2001, ApJ, 562, L15 
Tcplitz, H. I, Malkan, M. A., & McLean, I. S. 1998, ApJ, 506, 519 
Villar-Martin, M., Tadhunter, C, & Clark, N. 1997, A& A, 323, 21 
White, S. D. M., & Frenk, C. 1991, ApJ, 379, 52 
Windhorst, R. A., Keel, W. C, & Pascarelle, S. M. 1998, ApJ, 494, 

L27 



